KEK Preprint 95-14 
KEK-TH-436 
KEK-CP-27 



Finite Size Effects for the Ising Model Coupled to 

2-D Random Surfaces 



N. D. Hari DASsQ a , B. E. HanlonQ and T. Yukawa^] 

KEK Theory Group, 
1-1 Oho, 
Tsukuba, 
Ibaraki 305, 

Japan 



Abstract 

Finite size effects for the Ising Model coupled to two dimensional random sur- 
faces are studied by exploiting the exact results from the 2-matrix models. The 
fixed area partition function is numerically calculated with arbitrary precision by 
developing an efficient algorithm for recursively solving the quintic equations so en- 
countered. An analytic method for studying finite size effects is developed based on 
the behaviour of the free energy near its singular points. The generic form of finite 
size corrections so obtained are seen to be quite different from the phenomenolog- 
ical parameterisations used in the literature. The method of singularities is also 
applied to study the magnetic susceptibility. A brief discussion is presented on the 
implications of these results to the problem of a reliable determination of string 
susceptibility from numerical simulations. 
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1 Introduction 



Numerical simulations are important for the study of nonperturbative effects which are 
difficult to handle analytically. In such simulations one necessarily has to work with 
systems of finite size. The effects of finite size manifest themselves as systematic errors 
in measurements. These corrections are generically difficult to estimate as they involve 
dynamical details. A practical approach is thus often resorted to by applying phenomeno- 
logical parametrisations to finite size effects pj. 

In this paper we address the question of finite size corrections by considering the 
model of Ising spins coupled to two dimensional random surfaces. Since the Ising spin 
case, with central charge cm = 1/2, is exactly solvable by the method of matrix models[|], 
H U it provides an important test case for the efficacy of numerical simulations. An 
understanding of the nature of finite size corrections in this case may be useful in probing 
the unknown region beyond cm = 1. 

We approach this issue numerically and analytically by exploiting the parametric so- 
lution of the two matrix models^, [|]. In the numerical approach the inherent quintic 
equations are solved recursively from which the fixed area partition sums are extracted 
for various spin couplings. We also extend this approach to the case of nonzero magnetic 
field and determine the magnetic susceptibility scaling laws. The numerical analysis is 
augmented by an analytic analysis of the free energy about its singular points. From 
this we suggest a general ansatz for the form of finite size corrections. We apply our 
results to different estimates of the string susceptibility and to the minbu technique 0, 
demonstrating the difference with the phenomenological estimates previousy employed]!]]. 
The method of singularities is also applied to the problem of the finite size behaviour of 
the magnetic susceptibility. It is found that these agree with the usual finite size scaling 
laws0. 
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2 Numerical Solution of the Two Matrix Model 



Matrix models are solved by the method of orthogonal polynomials whereby a parametric 
solution for the free energy is obtained. The physical couplings are related to the param- 
eter so introduced by what we shall call a "constraint" . From such a solution the relevant 
critical exponents can be derived. We will restrict our analysis to surfaces topologically 
equivalent to S 2 . 

For pure gravity (the one matrix model) the constraint equations are quadratic and 
cubic for the case of quartic and cubic interactions respectively. The constraints are thus 
solvable by radicals allowing the constraint to be inverted. From this it is possible to 
derive a closed form for the free energy in the coupling g. For a quartic interaction the 
series solution for the free energy is given by[[|: 

(-12y)"(2n-l)! 
^ n\(n + 2)\ ' 1 1 

which can be written in the form T = J2n Z n g n , Z n being the fixed area partition sum for 
random surfaces. By exploiting the asymptotic expansion for the gamma function it is a 
simple matter to explicitly derive the asymptotic form of the fixed area partition sum: 

(-1)"(48)" , 25 

Zn 7^ ( ^ } • (2) 



The leading order behaviour corresponds to that originally predicted by KPZ|| and has 
the generic form Z n ~ exp(/in)n~ b , where the string susceptibility is defined as 7 = —6+3, 
and \x is the cosmological constant. Similar results can be obtained in the case of cubic 
interactions as well as for models with tadpole and/or self energy contributions removed 0]. 

To represent models in which Ising spins have been coupled to the, discretised, random 
surface a two matrix model is required. The method of solution follows as before. Again 
the free energy is given as a parametric solution in terms of a parameter z. However, we 
now have the difficulty of a constraint that is a quintic0, which for quartic interactions 
is given by ( where c = exp(— 2(3) and (3 = l/T): 



The general quintic is not solvable by radicals. As such a simple closed form for the free 
energy cannot be derived, as was possible in the case of pure gravity, and consequently 
the asymptotic form of the fixed area partition sums cannot be easily extracted. 

In this section we will describe our numerical method of solution to this problem where 
we seek a power series expansion for the free energy in terms of g. We will describe in 
detail the case for quartic interactions. The solution in the quartic case contains all the 
essential elements, while the cubic case[Q is both algorithmically and algebraically more 
cumbersome. 

To avoid the logarithm in the solution for the free energy we study dT jdg which is 
equivalent in this instance to dT jdg yielding: 

dT r 1 1 f z dt , . 1 f z dt . 
ag 2g g z Jo t g l Jo t 

We wish to represent this as a simple polynomial in z. To remove the inverse powers of 
(1 — 3z) which arise we employ the quintic constraint (|3]). dT jdg is then a polynomial in 
z of degree 8 with coefficients depending on the physical parameters c and g. This can 
further be reduced to a polynomial of degree 5 by using the constraint to express z 8 ,z 7 
and z 6 as quintic polynomials in z. 

To solve the quintic @ numerically we develop an efficient algorithm to handle the 
multiple sums inherent in a power series solution of this expression. We require an algo- 
rithm which avoids recalculation. We begin by explicitly expressing each power of z as a 
power series in g: 

z = $>i(n) g n ] z 2 = J2 oa(n) 2 5 = ^> 5 (n) # n . (5) 

Using (H), and the fact that 02(1) = 03(1) = 04(1) = 05(1) = it is a simple matter to 
extract ai(l) = 1/(1 — c 2 ). Indeed, a recursion relation for a\(n) for n > 2 can be derived 
from (^) which has the general form: 

oi(w) = 61 a i( n ~ !) + 6 a 2 (n) + £ 3 a 2 {n - 1) + £4 03(71) + £5 04(71) + £ 6 o B (n) . (6) 

where the £'s are the coefficients appearing in the quintic. The £12(71) 's, as(n)'s etc can 
ultimately be described in terms of the ai(n)'s. For each power of z the coefficients are 

4 



dependent on a calculation of those of lower orders such that for example 



N-l 



a 2 (N)=Y,a 1 (n)-a 1 (N-n) . (7) 
1 

Thus in order to calculate a 2 (N) we need to find ai(k) up to k = N — 1. Likewise, 
a 3 (N) = J^i' 1 o,i(n) ■ a 2 (N — n) implies that a calculation of a 3 (N) requires a 2 (N — 1), 
which requires a\(N — 2) etc. 

By this approach the calculation time is significantly reduced ( in this case the cal- 
culation time goes as ~ N 2 as opposed to ~ iV 6 for a naive approach to multiple sums). 
Substituting into (|J) and integrating with respect to g we can recover the fixed area 
partition sums up to large orders in a convenient time. 

We must, however, be careful to consider the exponential growth in the fixed area 
partition sums due to the cosmological term. To overcome this problem we estimate 
the cosmological term by studying the fixed area partition sums up to areas allowed by 
machine limits. The fixed area sums are then scaled by scaling all the a(N) by this 
estimate, so that 

a (iV) - e -» N a(N) ; 6 - e"^i , 6 - e""6 . (8) 

We found that in this way sufficiently accurate estimates of the cosmolgical constant could 
be obtained for relatively small values of n(~ 200) which could then be used to extract 
the scaled fixed area partition sums for even large n ~ 100000 easily. This makes it 
possible to investigate regions corresponding to those typical of numerical simulations via 
dynamical triangulation, and far beyond. 

To investigate the magnetic susceptibility we follow the same general prescription as 
above. With the introduction of a magnetic field two different coupling constants arise. 
In a perturbative solution to lowest order in H the constraint (|3|) becomes 

z z 2 H 2 

q(z,H) = - — - c 2 z + 3c 2 z 3 + — — . (9) 

yV ' ; (l-3z) 2 (1 -32) 2 (1 + 3z) 2 V ; 

Since we require solutions in the limit H — > the power series expansions for z, z 2 , ... etc 
in g remain defined as before by (^). The high and low temperature phases corresponding 



to the singularities of the free energy with respect to g are defined by the expression 
g'(z, H = 0) = 0. This has five solutions of which only two are physical for < c < 1: 

Zo = —1/3 (Low temp, phase) and z = Zq(c) (High temp, phase) (10) 

The critical temperature is at c = 1/4. 

The magnetic susceptibility is essentially given by the second derivative of T with 
respect to H at H = 0. In particular we find 

We can express this as before as a power series in g. We thus have an expansion of the 
form K(z) =Y1 K n g n . The fixed area magnetic susceptibility is then given by 

(12) 

This is valid at or above the critical temperature. Below the critical temperature we must 
account for the spontaneous ordering of spins so that the fixed area magnetic susceptibilty 
is given by 

Xn = ^-n<a> 2 , (13) 
nzj n 

where < a > is the spontaneous magnetization. 



3 Results from the Numerical Analysis 

We first exhibit the finite size effects in string susceptibility. For the case at hand 7 
is known exactly: 7 = —1/3 at the critical temperature and 7 = —1/2 off the critical 
temperature. We estimate the string susceptibility at finite area, 7 es t, by a suitable ratio 
of fixed area partition sums designed to cancel the cosmo logical constant: 

7est = {ln(^^)/ln(l-l)} + 3 . (14) 

This is a "local" estimate of 7 in that it involves neighbouring partition sums. Clearly, 
however, there are many different ways in which to extract such an estimate. 
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We present the results of 7 cs t for quartic interactions at the critical temperature in 
Fig(0) and simply note that the cubic case demonstrates the same behaviour. Indeed we 
observe the following important common features: 

(i) For large areas 7 es t approaches the theoretical expectation both on and off the critical 
temperature. 

(ii) Off the critical temperature 7 est approaches —1/2 rapidly from above. At the critical 
temperature 7 est drops below— 1/3 and then slowly converges towards the theoretical value 
from below. 

(iii) For the region at the critical temperature where 7 est is less than —1/3 the effects of 
finite size are greatest in the range approximately bounded by 200 < n < 2000. 

That a qualitative difference on and off the critical temperature should appear is 
consistent with the expectation that finite size effects will be influenced by large Ising spin 
correlations at the critical point. However, such a difference has not previously been taken 
into account. We observe that this is, in fact, an important factor which cannot be ignored. 
Significantly, the effects of such finite size corrections appear most pronounced in regions 
where previous numerical simulations have concentrated their estimates of 7, which most 
typically deal with simulations of the size ranging around n = 1000 ~ 2000 [I]. We see 
that finite size effects at the critical temperature do not simply decrease with increasing 
size but exhibit an important nonlinearity with greatest effect around the simulation sizes 
chosen. In addition we note that, since a range of values for 7 es t are recovered, it follows 
that we can engineer to find many values of 7. Indeed, at the critical temperature 7 es t 
actually passes through the correct value even for small systems. It appears then that 
extraction of 7 from small area studies can be misleading. 

As we pointed out above, these results are exhibited for a particular choice of form for 
the estimation of 7. However, as will be demonstrated below, while some of the features 
may change, our observations in general remain sound. 

The magnetic susceptibility results reproduce those expected from Liouville theory fll0| 
as well as from standard scaling analysis ||. In particular we find that 

Xn — > constant (High temp, phase) 
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n 2 ' 3 (on the critical temp.) , 



(15) 



while in the low temperature phase 

Xn + n < o > 2 —> n . (16) 

4 Singularity Analysis of Finite Size Effects 

We wish now to analytically find the asymptotic form of the fixed area partition sums. 
This is made possible by the observation that the large n behaviour of Z n is dominated 
by the singular points of the free energy. As before we will deal explicitly with the case 
of quartic interactions. We motivate the analysis by applying this approach to the simple 
case of pure gravity. 

The parametric solution for the one matrix model with quartic interactions is given 
by! 

^=-ilns+J-(*-l)(9-s) with g( z )^ l —± = g . (17) 

The singularity of T with respect to g is determined by the condition g'(z Q ) = for which 
zq = 2, so that g c = —1/48. The constraint in (|PT|) can easily be inverted to yield 

(z - z ) = a{g - g c ) 1/2 + b(g - g c ) + (18) 

where for brevity we have not displayed the values of a, b, .... This in turn can be used to 
generate an expansion for T about the critical point g c . Being careful to retain sufficient 
terms in the expansion the nonregular terms are found to be 

. , 12283^3, x , /2 1769472^, N7/2 

Hg) - 9c?' 2 + 7r-^{g - 9c) 7/2 + 19 

5 7 

The required asymptotic form can be obtained by employing the binomial expansion 
for (g — g c ) a arid subsequently the asymptotic expansion for the gamma function. Writing 
(9 ~ 9c) a = En 9 nz n , one finds 

Z% ~ -(-^) a ~"(-l)" Sm(7ra) r(a + l)n- {1+a) exp(a(l + a)/2n) + (20) 

7T 
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The term with a = 5/2 is seen to reproduce the leading order behaviour found in @. 
We can also find the leading order corrections which involve contributions from both the 
(g — g c ) 5 ^ 2 and (g — g c ) 7 ^ 2 terms: 

7 F ieadin gn ,35 1769472 5 1 ^(9/2)1 

Zn ~ Zn (1 + 8^ ( ~ 7~ )( lW ( 48 } f(772)n + -- ) (21) 

which is exactly that expressed in (0). We thus have a method by which finite size 
corrections may be derived without relying on a closed form for T . 

We now apply this approach to the two matrix model. From the parametric represen- 
tation of the free energy it follows that T'{zq) = dJ- '/dz\ Zo = 0, where zq is given in (|T0|). 
The behaviour of the second derivatives with respect to z are: 

Off crit. temp. : ^"{zq) ^ 0, g"{z ) ^ , 

On crit. temp. : T"{z Q ) = 0, g"{z ) = , (22) 

where the non-vanishing of jF"(z ) off criticality is crucial to obtaining the correct scaling 
laws. Hence we differ from[]3| in this respect. We must thus consider two regions, 
(i) The case off the critical temperature, c ^ c cr ; t : As in the case of pure gravity we 
invert the constraint (^) to give an expansion for (z — ^o) in (g — g c ): 

(z - z ) = a{g - g c f/ 2 + b(g-g c ) + (23) 

Taylor expanding T around z and substituting (|23| ) generates a series of both regular 
and nonregular terms in (g — g c ). Formally, the lowest exponent of the nonregular terms 
is 3/2 but this term vanishes owing to the relation 

F'(z )ab+^F"\z ) = , (24) 
o 

explicitly requiring that T"(zq) ^ 0. The nonregular terms contributing to JF are thus 

T{g) ~ Al{g - g c f 2 + A2(g - g c ) 7 l 2 + .... (25) 

so that 

Z n - (-g c )- n n- 7/2 {l + — + ™ + •••} , (26) 

n n l 
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again exhibiting the same basic form as that for pure gravity. Some representative values 
for Bl with quartic interaction are 

Bl = -72.69 at c = 0.20 ; Bl = -8.76 at c = 0.36. (27) 

(ii) At the critical temperature c = c cr i t : Since now g"(zo) = the expansion of (|3|) 

takes the form (g — g c ) — (z — z ) 3 g'" (z ) / 6 + (z — zo) 4: g IV (zo)/24 + , which after 

inversion gives : 

(z - z ) = a(g - g c f' z + b(g - g c ) 2 ' 3 + d(g - g c ) + (28) 

Furthermore, since now T"{zq) = 0, the Taylor expansion for T around zq starts at 

(z — -2 ) 3 . Substituting ( p8|) into this expansion for T the coefficients of powers of (g — g c ) 

conspire so that the (g — g c ) 4 ^ 3 and (g — g c ) 5 ^ 3 terms are absent. The leading singular 

behaviour for T is thus given by T ~ (g — g c ) 7 ^ 3 . Consequently, following the same steps 

as before, the leading behaviour of Z n will take the form Z n ~ (— g c )~ n n~ 10 / 3 , from which 

it follows that the string susceptibility is given by 7 = —1/3. The corrections can be 

similarly calculated for which we quote the results: 

• • • ~ ^win^r, 0.4287 3.08 1.2980 
Quartic interaction : Z n = Z„ mg {l + — + ...} 

^ Pa H in ,r 0.286 3.05 0.936 , . , 

Cubic interaction : Z n = ^ eadmg {l + ttt + •••} • 29 

The crucial observation here is that the next to leading order correction goes as 1/n 1 / 3 
rather than 1/n which is the case with pure gravity. We have verified that these corrections 
reproduce the observed finite size effect from the numerical analysis (see Fig(|2|)). We have 
thus isolated the fundamental difference between these two cases. It is now clear that 
assuming 1/n type corrections both off and on the critical temperature is not justified. 

For the magnetic susceptibility the important expression was given in (|TT| ) which we 
can express as the sum of two integrals K(z) = I\ + I<i- We see from (|^) that at the 
critical point z the expression g(z,H) is singular. Consequently both I x and J 2 are 
singular. However, the sum I\ + 1 2 is regular at z . 

As with the free energy we can extract the critical behaviour of the magnetic suscepti- 
bility by taking account of its behaviour about the singular point. Expanding K(z) about 
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zq we find that, in the high temperature phase as well as at the critical temperature, 
K(z) ~ at(z — zo) 2 + (3(z — ^o) 3 + ...., while in the low-temperature phase, K(z) ~ (z — Zq). 
We know from (|23| ) that off the critical temperature {z — Zq) ~ a(g — gc) 1 ^ 2 + b(g — g c ) + .... 
Consequently, the leading nonregular term is K(g) ~ (g — g c ) 3 ^ 2 , where we have verified 
that this term does not vanish as was the case with the free energy for which we know 
that T ~ (g — g c ) 5 ^ 2 - Expanding in g n it follows that 

K n 

Xn = — ~ ► constant . (30) 

Similarly at the critical temperature we have [z — Zq) ~ a(g — gc) 1 ' 3 + ••• so that K(g) ~ 
(g — g c ) 2 ^ 3 - We know that T ~ (g — g c ) 7 ^ 3 so that we find 

Xn ~ n 2 ' 3 . (31) 

In the low temperature phase, it follows that K(g) ~ (g — g c Y^ 3 - Consequently, 

K„ 



nZ r , 



n . (32) 



We can thus account for the behaviour expected from standard Fisher-scaling theory || 
as well as from Liouville theory flQ|. 



5 Alternate Estimates of 7 



We now investigate alternate definitions for the estimation of 7. From Z n one can intro- 
duce an obvious such alternative: 

ln(l + finite size corrections) 

7alt. est — 7exact H , / n > (33) 

m(n) 

from which it follows that 

Cl 



Quit • 7alt. est ~ 7exact T 



ln^n 1 / 3 



C ^ C C nt '■ 7alt. est ~ 7exact + , , \ , (34) 

m(n)n 
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where c\ > and a\ < 0. According to this definition the finite size behaviour predicted is 
opposite that found from the numerical results for 7 est . That is, off the critical temperature 
7ait. est approaches —1/2 from below while on the critical temperature it approaches —1/3 
from above. There is no real inconsistency here as there are many ways in which to 
estimate 7, and the different estimates are only required to coincide asymptotically. For 
example, Brezin and Hikami|| use an estimate based on Pade approximation where 

_ o n(l + n)(f n - / w _i) 

7Pade est — & 77—. w 7 , (.OO) 

(l + n)f n - n/ n _i 

and /„ is a Pade approximant to the ratio Z n jZ n -\. We observe that each choice has a 
different finite size behaviour. We can explicitly demonstrate this where in addition to 
those results in ([34]) we have that on the critical temperature, where Z n ~ i^ eadmg (l + 
Ci/n 1 / 3 + c 2 /n + ...) and C\ > 0, 

1 4 Cl 5c? 24/3 + c 2 

7est ~ "3-9^ + 9^ ~ + ^ 

1 4ci 5c? -2c 2 /3 - 2c 2 - 36 + 6 2 

7Padeest ~ "3-9^+9^ + n + "" (36) 

where b = 3—7, while off the critical temperature where Z n ~ Z^ admg (l+ai/n+a2/n 2 +...) 
and ai < 0, 



1 3a? — 6a 2 
7est ~ ---2ai/nH h .... and 

2 n z 



1 20! + 3o - o 2 

7Padeest ~ ~ ~ H 

2 n 

where we note that — 2ai — 36 + 6 2 is positive and that 7 es t is in fact analytically consistent 
with the numerical results. The situation is thus similar to the scheme dependence in 
renormalisation. Any particular prescription for 7 is suitable but comparison of 
different definitions is not a meaningful exercise. Thus different models should be 
compared using the same definition for 7 est . From these analyses we can summarise the 
general structure of finite size corrections pertinent to any choice of estimate: 

(i) They are not always parameterisable as 1 + a/n + . . . . 

(ii) The parameterisation is dependent on the expected value of 7. 
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(iii) If 7ex P ectcd = q/p for q,p G Z then in general the finite size corrections will be 
parameterised as 1 + a/n 1 ^ + .... This is our general ansatz for the form of finite size 
corrections. 

(iv) Some of the coefficients in the expansion may be zero, for instance in the Ising case 
at the critical temperature we have the expansion 1 + ci/n 1 ^ 3 + C2/n + c^/n 4 ^ 3 + c^/n 2 + ... 



6 Some Comments on Minbu Analysis 

A particular approach to extracting the string susceptibilty from numerical simulations is 
given by measuring the distribution of minimum neck baby universes (minbus) on random 
surfaces of a given fixed total area A. An estimation for the average number of minbus 
of size B, ua{B), was given by Jain and Mathur || and takes the form 

For quartic interactions we should substitute 4 for 3. It follows that (to leading order) 

\n(n A (B)) = constant + (7 - 2) ln(5(l - — )) . (39) 

Thus by measuring the slope of this function we can numerically determine 70. 

The form of finite size corrections for pure gravity being (1 + a/n+ ...) it follows that 



the leading finite size correction to (|39|) is 



ln(n A (S))=constant + (7-2)ln(S(l-^)) + (7-2+a)ln(l+— -L— —) + ... , (40) 

We know from @ that for pure gravity with quartic interactions the finite size parameter 
has the value a = —25/8. Plotting the minbu distribution ( [40"1) and the exact minbu 
distribution calculated using (1) in Fig@) we see that this correction gives excellent 
agreement with the exact result. Applying this approach to the Ising case off the critical 



temperature, where the finite size corrections were given in (|27|) for the high and low 
temperature phases, we find that these finite size corrections do not accurately mimic the 
exact minbu distributions, Fig(Si). Clearly, higher order corrections are more important 
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in the Ising case. This behaviour is again evident at the critical temperature, where we 
must also account for 1/n 1 / 3 type corrections. Again we find that retaining higher order 
corrections is necessary in order to obtain a reasonable fit. 

Since minbus are measured over a range of volumes finite size effects are unavoidable 
for small minbus. As these finite size effects are hard to extract for cases where either 7 
is not known or where 7 = q/p with large integer p, it appears that applying the minbu 
technique to extract meaningful estimates for 7 from simulations in these interesting cases 
is fraught with difficulties. 

7 Conlusion 

A simple parameterisation of finite size effects is a natural first approach to analysing 
numerical simulations. As we have shown, however, the actual parameterisation is non- 
trivially dependent on the exact value of the string susceptibility. On this basis we have 
proposed a general ansatz for the form of finite size corrections. A possible algorithm 
then is to make a best guess for 7 and to fit this to the observed data with the finite size 
corrections corrrectly included. By performing a x 2 analysis the best fit for a particular 

7 could be recursively searched. 

We have compared different approaches to estimating 7 and demonstrated that at- 
tempts to compare estimates from different definitions can be misleading. These consid- 
erations become relevant if we wish to extract reliable numerical estimates beyond the 
Cm = 1 barrier where it is known that large logarithmic corrections also play a role. A 
clear understanding of the functional form for the finite area estimates are indispensable 
in these cases. 
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Figure 1: 7 es t at the critical temperature for quartic interactions. 
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Figure 2: Comparison of the analytic and numerical estimates of 7 at the critical tem- 
perature with quartic interactions. The analytic graph includes all the corrections given 
in (29). 
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Figure 3: Numerically and analytically derived minbu plots for (i) pure gravity and (ii) 
the Ising case off the critical temp, (c = 0.36). The calculated values of \n(nA(B)) in (ii) 
have been shifted by a constant for clarity. All are with quartic interactions. 
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